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We present a study of the prospects for coherence preservation in solid-state spin qubits using 
dynamical decoupling protocols. Recent experiments have provided the first demonstrations of mul- 
tipulse dynamical decoupling sequences in this qubit system, but quantitative analyses of potential 
coherence improvements have been hampered by a lack of concrete knowledge of the relevant noise 
processes. We present simulations of qubit coherence under the application of arbitrary dynam- 
ical decoupling pulse sequences based on an experimentally validated semiclassical model. This 
phenomenological approach bundles the details of underlying noise processes into a single exper- 
imentally relevant noise power spectral density. Our results show that the dominant features of 
experimental measurements in a two-electron singlet-triplet spin qubit can be replicated using a 
l/cj^ noise power spectrum associated with nuclear-spin- flips in the host material. Beginning with 
this validation we address the effects of nuclear programming, high-frequency nuclear- spin dynamics, 
and other high-frequency classical noise sources, with conjectures supported by physical arguments 
and microscopic calculations where relevant. Our results provide expected performance bounds and 
identify diagnostic metrics that can be measured experimentally in order to better elucidate the 
underlying nuclear spin dynamics. 



I. INTRODUCTION 

Spin qubits in semiconductor materials are a leading 
experimental realization of a controllable, scalable quan- 
tum system in the solid-state^. These devices hold nu- 
merous advantages relative to other qubit implementa- 
tions; First, the use of a semiconductor platform pro- 
vides benefits in terms of scaling - leveraging decades 
of development in advanced micro- and nanofabrication, 
large-scale semiconductor-based systems appear techni- 
cally feasible. Second, experimental demonstrations of 
coherent control of semiconductor spin qubits allow op- 
erations in the nanosecond and sub-nanosecond range^ ^, 
ultimately permitting rapid logic operations at the phys- 
ical qubit level. Finally, measured electron-spin relax- 
ation times (Ti) in semiconductor nanostructures exceed 
several seconds^, placing the limit on coherence time 
^rj.{max) ^ up to 10 orders of magnitude larger than 
demonstrated gate times, an important DiVincenzo cri- 
teriorP. 

Many realizations of semiconductor spin qubits have 
been proposed and demonstrated experiment all}'^. Ad- 
ditionally, experimental control fidelity has progressed to 
the point where probing studies of the coherence limits of 
spins in semiconductors are possible^. We focus on the 
two-electron Singlet-Triplet (S/T) spin qubit in GaAs^ 
as it represents one of the most advanced solid-state spin 
qubits available. This system employs the symmetric and 
antisymmetric spin configurations of a pair of spin- 1/2 
electrons to form a qubit basis. Electron pairs are con- 
fined in a lithographically pattered nanostructure based 
on a two-dimensional electron gas and confining electro- 
static gates. Decades of research in mesoscopic physics 



now permit experimentalists to isolate electron pairs, en- 
act controlled exchange-based interactions, and perform 
a strong projective measurement via spin-to-charge con- 
version and use of a proximal charge detectoi^. 

The strengths of the S/T qubit, however, come at sig- 
nificant cost in terms of decoherence processes. Most sig- 
nificantly, the localization of the qubit within a nanopat- 
terned GaAs-heterostructure host material introduces 
noise sources due to fluctuating nuclear spins^^ and lo- 
calized charge centers^^. A series of microscopic theoreti- 
cal studies and experiments have shown that in these de- 
vices the effects of nuclear spin dynamics dominate mea- 
sured qubit coherence. To date, typical free-induction- 
decay (FID) T2^^^^ times have been measured ~20-50 
j^^ 2 | i2 | i3 |^ and spin echo experiments (with exchange off 
have yielded coherence times j"!^^^^^^ ~ 30 ^s. Recent 

experiments have extended plf'^^^ to hundreds of ns us- 
ing nuclear programming^^, and have demonstrated long 
coherenc e time s using multipulse dynamical decoupling 
sequenceJI^'I^ to mitigate the effects of dephasing. At 
this time, however, an efficient and practical model pro- 
viding quantitative insight into the ultimate prospects for 
the suppression of decoherence-induced error accumula- 
tion remains elusive. 

In this manuscript we use an experimentally validated 
theoretical model for error accumulation under the ap- 
plication of dynamical decoupling pulse sequences to ac- 
complish four main goals: 1) Provide a general phe- 
nomenological framework for spin-qubit error accumu- 
lation in terms of noise power spectral densities by com- 
parison with experimental measurements, 2) Provide in- 
sight into the spectral characteristics of dominant noise 
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processes, 3) Understand the prospects for and hmita- 
tions of dynamical decoupling in extending the coher- 
ence of spin qubits, and 4) Identify the most promis- 
ing technical approaches to improving qubit coherence. 
Our numerical simulations first suggest the coherence 
times achievable under dephasing due to fluctuating nu- 
clear spins. Calibrations of the noise power spectral den- 
sity, 5/3 (cj), for angular frequency cj, are performed us- 
ing measurements of T^f^^^^ and the relative coherence 
time for an n-pulse CPMG cycle relative to the FID 
time, j"^^^^^) Our calculations allow com- 
parisons of the scaling of measured coherence times with 
pulse number, n, and studies of the influence of differ- 
ent spectral components of S(3{uj)^ particularly the form 
of the high-frequency cutoff. Further, we are able to 
study predicted error rates under dynamical decoupling 
in the high-fidelity regime, vital for quantum computa- 
tion. This model allows detailed comparison with exper- 
iment and provides predictive power for expected qubit 
coherence times in the presence of realistic noise sources. 

The remainder of this manuscript is organized as fol- 
lows: Section [Til introduces the theoretical model used for 
the simulations and describes the relevant experiments. 
This is followed by detailed simulations in Section |III| 
of predicted coherence in the presence of a noise 
power spectral density. More complex power spectra and 
the role of a high-frequency noise cutoff are introduced 
in Section IIVI The effects of noninstantaneous and im- 
perfect control pulses are described in Section [V| and the 



manuscript concludes with discussion in Section VI 



II. 



DEPHASING NOISE IN A DYNAMICAL 
DECOUPLING FRAMEWORK 



We consider a simple theoretical model based on pre- 
vious work^^ and validated by experiment^^ in or- 
der to describe dephasing in a generic qubit system. We 
address only phase randomization due to classical envi- 
ronmental noise, described by a Hamiltonian written as 



F= i[!^ + /3(t)]az, 



(1) 



where Vt is the unperturbed qubit splitting, /3 is a time- 
dependent classical random variable, and cfz is a Pauli 
operator. For our purposes a full quantum mechanical 
treatment is not required, although it is possible follow- 
ing previous work^^^. Further, as spin lifetimes in GaAs 
can exceed 1 s, current experiments are dominated by the 
effects of dephasing. 

The term I3(t) captures environmental fluctuations 
that produce an effect ive, fluctuating magnetic field on 
the qubit systeni^^^ all qubit-specific features of the 
model are captured through This term may appear 
due to external fluctuating magnetic fields, or intrinsic 
processes such as nuclear spin flips. 

In the frame rotating at Vt^ j3(t) produces a random 



phase shift between the qubit basis states that on aver- 
age leads to a ~ 1/e decay in coherence when the root- 
mean-squared phase accumulation is ~ V^tt. The char- 
acteristic timescale for this process is known as r^, and 
in the absence of relaxation is equivalent to T2, generi- 
cally known as the decoherence time. For the remainder 
of this manuscript we refer to the 1/e coherence time of 
the qubit as T2. 

From an experimentalist's perspective it is useful to 
characterize P{t) in the frequency domain, using the noise 
power spectrum Sp{uj)^ the Fourier transform of the two- 
time correlation function of 
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m+T)m)dT. 



(2) 



Here uj is angular frequency. The influence of the noise 
term after time r enters the measure of qubit coherence 
for a superposition state in the equatorial plane of the 
Bloch sphere by writing 

-X(r) 



W{T) = \{cJy){T)\ 



(3) 



where angled brackets indicate a quantum-mechanical ex- 
pectation value, the overline indicates an ensemble aver- 
age, and 



X(r) 



00 

11 



F{ujT)duj. 



(4) 



In the expression above, F{ujr) is known as a filter- 
functiod- ^^ l -'^^ I ^^ I ^^ J which encapsulates the experimental 
conditions under which qubit coherence is measured. 
Since the filter function enters the coherence integral as 
a multiplicative factor of ^'/^(a;), small values of F{ujr) 
where 5'/3(a;) is large will lead to small values of x{^)^ 
and hence coherence W{r) ~ 1. 



A. Free- Induct ion Decay 

In a Free-Induction-Decay (FID) experiment the filter 
function takes the form 



F{ujr) = 4sin^(cjr/2), 



(5) 



such that qubit coherence is set by the integral over all 
spectral components of 5'/3(cj), when the small angle ap- 
proximation is valid. This intuitively makes sense, as in a 
FID experiment the precise spectral characteristics of the 
noise are immaterial - only the net average phase accu- 
mulation matters, consistent with the fact that the phase 

accumulation in the time domain is given 

Low-frequency fluctuations enter into an ensemble av- 
erage measure of decoherence as they produce shot-to- 
shot quasistatic phase offsets at the conclusion of an ex- 
periment. When averaged over many experiments these 
phase offsets produce a decay in coherence. This form of 
decoherence is sometimes referred to as a measurement 
of T2* in a time-ensemble-averaged fashion. 
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B. Dynamical Decoupling Sequences 

Dynamical decoupling is a technique derived from the 
Nuclear Magnetic Resonance (NMR) community that 
has been proven useful for suppressing decoherence in 
a quantum informatic settin^^ ^-^ l ^^ * ^^" ^. In a dynam- 
ical decoupling framework, the free evolution of a qubit 
is broken into discrete time periods between which time- 
reversing operations are applied, effectively decoupling 
the qubit from its fluctuating environment. This ap- 
proach is based on Hahn's discover};^^ of the spin echo 
in NMR for the mitigation of inhomogeneous dephasing, 
but applies equally well to the case of a single-spin. 

The quantitative effect of dynamical decoupling is cap- 
tured by noting that the action of intermittent applica- 
tion of cFy operators ("tt" -pulses about the F-axis) leads 

the phase accumulation, e , to be broken into 

segments corresponding to the interpulse periods, with 
the sign of the phase accumulation alternating in succes- 
sive periods. Uhrig^^ and Cywinski et al.^^ showed that 
for any n-pulse sequence one may account for this by 
writing a time-domain filter function, ^n(^), with values 
±1, alternating between each interpulse free-precession 
period. We then write the filter function in the frequency 
domain, F{ujt) = \yn{ujr)\'^, where yn{ujr) is the Fourier 
transform of the time-domain filter function, and r is the 
total sequence length. Again, the filter function describes 
phase accumulation in the frequency domain under the 
application of a dynamical decoupling pulse sequence. 
For an arbitrary n-pulse sequence we may thus write 

F{uT) = \yn{^r)\^ 
= |1 + (-i)n+ie^^^ + 2 ^ (-Ipe^^^-^^P (6) 

where 5jT is the time of the j*^ ttx pulse. For conve- 
nience we assume instantaneous pulses, but we address 
this point in a later section of this manuscript (Sec- 
tion |V|). 

Using this model we may therefore determine the av- 
erage dephasing error expected for any qubit given a well 
defined S 13(00) and the analytically defined filter function 
for an arbitrary applied dynamical decoupling sequence. 
This is an extremely efficient and practical model al- 
lowing all details of system-bath interactions (including 
external as well well as intrinsic microscopic noise pro- 
cesses) to be represented through the form of S(3{uj). This 
construction has been quantitatively validated by exper- 
iments using trapped atomic which good 
agreement was found between experiment and theory for 
a wide range of S(3{uj) and applied pulse sequences. 



C. Dynamical Decoupling Sequences of Interest 

In this study we focus on the potential error- 
suppressing capabilities of two well studied sequences 



that have attracted significant attention in the last sev- 
eral years, CPMG and UDD. These sequences form an 
instructive set to study due to their relative ease of im- 
plementation and distinct spectral-filtering characteris- 

The n-pulse CPMG sequence is an extension of the 
spin-echo sequence and has been widely utilized in NMR 
and quantum information experiments^^ The se- 
quence has been shown effective at suppressing phase 
randomization when noise processes are dominated by 
low- frequency components (e.g. S 13(00) ex 1/oj^P^. An n- 
pulse CPMG sequence has all pulses evenly spaced, with 
the first and last free precession periods half as long as 
the interpulse free precession periods. For the purposes 
of this study we do not differentiate between the Carr- 
Purcell and CPMG variants of the multipulse spin echo 
(CPMG is more effective at suppressing pulse rotation er- 
rors and most efficiently suppresses dephasing only when 
the initial Bloch vector is aligned with the applied trans- 
verse field). 

Uhrig analytically derived^ an n-pulse sequence in 
which the first n derivatives of yn{ojr) vanish for ojt = 0, 
exclusively through manipulation of the relative pulse lo- 
cations, Sj. The sequence is constructed using tt pulse 

locations determined analytically as Sj^^ = sin^[7rj/(2n + 
2)], for an n-pulse sequence. The resulting sequence, 
UDD, has been shown to be extremely effective at sup- 
pression of high-frequency noise with a sharp spectral 
cutoff, appropriate for noise that may be present in, 
e.g. a spin-boson model. The performance of UDD 
has been studied extensively theoretically and experi- 
mentall3n™llEa. 

Both the CPMG and UDD sequences are capable of ex- 
tending the qubit coherence time, but UDD has proven 
especially efficient at suppressing dephasing-induced er- 
ror at times short compared to the qubit's 1/e decay time. 
This regime is especially important for quantum informa- 
tion processing where constraints imparted by quantum 
error correction necessitate exceedingly low single-qubit 
error rates. 

Additionally, the relative differences in the form of 
the filter functions of CPMG and UDD, and the com- 
mensurate performance differences under S 13(00) dom- 
inated by either low-frequency contributions or high- 
frequency spectral components, make the use of these two 
sequences an interesting diagnostic tool for noise spec- 
troscopy. Quantitatively characterizing the expected de- 
phasing error in a qubit under application of UDD and 
CPMG for different n and r will therefore assist in re- 
constructing the spectral characteristics of relevant noise 
sources. 



III. THE INFLUENCE OF l/u^ NOISE ON 
SINGLET-TRIPLET QUBIT COHERENCE 

A number of decoherence sources have been identified 
in S/T qubits, with dominant sources relating to the tem- 
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poral evolution of nuclear spins in the GaAs heterostruc- 
ture^ 10 40 ^ process known as spectral diffusion permits 
a zero-energy exchange process between distant nuclear 
spins, changing the Overhauser field experienced by the 
electrons, and thus causing dephasing. Statistical mea- 
surements of the S/T qubit's singlet-return-probability 
suggest the presence of an effective noise power spectral 
density 5/3 (cj) oc XjuP' . While these nuclei can be treated 
as a fully quantum mechanical bath, it is sufficient to 
treat the ffuctuating Overhauser field associated with the 
nuclei as a classical noise source producing dephasing. 

We therefore begin with a study of the infiuence of 
S^{uj) (X XjuJ^ on dynamical decoupling performance, 
informed by experimental observations. We bound our 
spectral range of interest uo / 2t: ^ [0.01, 10^] Hz, in order 
to prevent numerical divergences. The selection of these 
bounds is consistent with both experimental observations 
and general physical arguments. Typical experiments on 
S/T qubits show Tlf'^^^ ^ 20-50 ns, and maximum co- 
herence times with multipulse dynamical decoupling se- 
quences 2^^(^^^^i6) ^ 200/is. A noise contribution oscil- 
lating with frequency 0.01 Hz corresponds to a 100 s pe- 
riod, indicating that this contribution is static to roughly 
10~^ over even the longest experiments (~ 10 for a 
FID experiment). We have confirmed in the forthcoming 
simulations that reducing this low-frequency bound fur- 
ther negligibly impacts the results. 

Similarly the high-frequency cutoff is motivated by 
bounds of the so-called dynamical decoupling limit. If 
ffuctuations in the noise are sufficiently rapid compared 
to the smallest interpulse spacing in an n-pulse sequence, 
^max /27r ^ ^^in^ ^^eu the application of a multipulse 
sequence can only lead to an inc rease in net dephasing 
error. Experiments on S/T qubit^^^^ have used up to 
n ^ 20, indicating a value of Tmin ^ 5/is = (200 kHz)~^ 
for CPMG. The selected value of uOrnax/"^^ ^ ^min^ 
the form of S(3{uj) ensures that components at these fre- 
quencies do not contribute significantly. 



A. Calibrating the Noise Strength 

We assess the scaling factor, a, that serves to set the 
overall magnitude of the applied noise, writing S(3{uj) = 
auj~'^, such that for the minimum frequency defined 
within the spectral range above, S^{27r 0.01 Hz) = a. We 
determine the appropriate scaling by studying the calcu- 
lated FID signal for various values of a. Physical insight 
into the approximate values of a is derived by considering 
the phase evolution in a FID experiment. In the small 
angle approximation the ensemble averaged dephasing 
is given by W{t) = exp[—a{2/7r)t'^J^Si3{uj)duj] = 



exp 



(t/r: 



(F/D)y 



This corresponds to a Gaussian 



decay where all factors (other than t^) may be combined 
to give the 1/e decay constant, T2^^^\ 

The 1/e coherence time of the FID signal is reduced ex- 
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FIG. 1: Qubit coherence in the presence of l/cj^ noise, (a) 
Effect of noise scaling, a, on predicted T2^^^\ A value of 
a ^ 5 X 10 most closely reproduces the observed T2 
The dashed line represents the 1/e error rate, (b) Expected 
1/e coherence time for CPMG and UDD pulse sequences as 
a function of pulse number for different values of a. 



ponentially with a, as demonstrated in Fig. [T^. Here we 
show the accumulated error (dephasing) signal as a func- 
tion of the FID experiment time, for various values of a. 
Our measure of coherence is presented as (1 — W{t))/2 
a probability of qubit measurement in a selected basis 
state at the conclusion of the experiment. A value of 0.5 
corresponds to 50% probability of the qubit being mea- 
sured in \S) or |T) — complete dephasing. 

Comparison of simulations to experimental measure- 
ments gives an estimate of the noise scaling factor, a. 
First, we look to measured values of T2^^^\ and find 

values of a = 10^^ - 10^^ reproduce T!f^^^ 20 - 50 
ns with an approximately Gaussian decay. This range of 
a gives a predicted value of j^j^^^^^^ ^ 20 — 60 /is, corre- 
sponding to an improvement in coherence time of three 
orders of magnitude using only a single pulse, and con- 
sistent with experiments^ 

Further information about the appropriate selection of 
a can be derived by examining the calculated values of 
the 1/e coherence time for n-pulse dynamical decoupling 
sequences (^j^^^^^^) d.ndT^^^'^^) {Y'lg.^p). Compar- 
ison against experimentally measured values of T2^^^\ 

^iEcHo)^ ^iCPMG,)^ ^iCPMG.o)^ ^(CPMG,,) ^ ^^^^^^ 

our approximation to a ~ 5 x 10^^. We note explicitly 
that this assumption is limited by measurement impreci- 
sion and variability between experiments, as well as lim- 
itations in our model. 

The value of a extracted from our simple phenomeno- 
logical model can be compared with experimental mea- 
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surements of the fluctuating nuclear fleld believed to be 
responsible for the XjuJ^ form of the noise power spec- 
tral density. Direct measurements of the noise spectrum 
around 1 Hz, where it already shows l/cj^ behavior in- 
dependently leads to a « 2.2 x 10^^, within a factor of 
order unity of the value extracted from our simulations. 



B. Coherence with Multipulse DD 

The data presented in Figure [TJd provide information 
on the scaling of coherence time with n, and may be 
compared against full analytical calculations for spin- 
dephasing and experimental measurements. Qualita- 
tively, we see that the spin-echo, CPMG, and UDD se- 
quences all efficiently eliminate low-frequency contribu- 
tions to dephasing. Thus it is observed that application 
of even a single spin-echo pulse will signiflcantly extend 
qubit coherence^^, and the application of sequences with 
increasing values of n will provide diminishing returns. 
Data for flxed a are well approximated by a power law 



scaling 



(CPMGr,) 



(X 



or TJ 



(X 
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(indicated by an approximate linear scaling on a log-log 
plot). Best flt values of for the two sequences are 
^{CPMG) ^ Q QQY and = 0.620. These results 

are commensurate with analytical predictions of Ref. [19] 
in which = § was predicted in the presence of Gaussian 
noise. 

In these calculations UDD is shown to generally yield 
shorter values of the measured T2 relative to CPMG for 
this S(3{uj). This is consistent with experimental mea- 
surements in a variety of systems, and is expected from 
a simple examination of the fllter functions for CPMG 
and UDD; the UDD sequence trades an extremely small 
value of F{uot) at low uo for a slightly lower "turn-on" fre- 
quency, above which the fllter function does not suppress 
dephasin^^ 



IV. THE HIGH-FREQUENCY CUTOFF 

The results presented thus far, and the similarity be- 
tween these numerical calculations and data obtained 
by both experiments and more careful analytical ap- 
proaches lend weight to the utility of this phenomeno- 
logical approach. In this section we will examine the 
potential impact of more carefully modeling the high- 
frequency nuclear spin dynamics. The error accumula- 
tion in a FID experiment is again set by the integrated 
noise power spectral density, meaning that modiflcation 
of the high-frequency contributions to Sjsiuo) yields negli- 
gible changes. Unlike the discussion of the low-frequency 
cutoff above, it is anticipated that the specific details of 
high-frequency noise will significantly impact the perfor- 
mance of the spin echo and multipulse dynamical decou- 
pling 

We consider modification of S^iuj) to account for the 
high-frequency dynamics of the nuclear spin bath. Above 
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FIG. 2: The effects of a high-frequency cutoff on calculated 
spin echo signal, (a) Schematic of S{uj) including an high- 
frequency cutoff to account for the dynamics of nuclear flip- 
flops. The high-frequency cutoff, ujc is variable, as is the form 
of the rolloff above uJc • (b) 7 as a function of uJc for different 
forms of the high-frequency rolloff. (Inset) j'^^^^^^ a func- 
tion of uc for the same parameters, (c) Error probability as 
a function of free-precession time in a spin echo experiment. 
Data approximated with a best fit superexponential: No cut- 
off, 7=3; 12 dB/octave, 7 = 3.62; 24 dB/octave, 7 = 3.81; 
36 dB/octave, 7 = 3.84; Exponential, 7 = 3.65; Gaussian 



7 = 3.81. For these rolloff s 
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for fits using 7 = 4: 12 dB/octave, x -0.010; 24 dB/octave, 
X^=0.002; 36 dB/octave, x^=0.002; Exponential, x^=0.009; 
Gaussian x^=0.001. a = 5 x 10^^ 
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a correlation time rc = ^njujci it has been calculated 
that spectral diffusion processes are suppressed. This 
corresponds to a high-frequency cutoff, cjc above which 
S(3{uj) is reduced below scaling. We study the effects 
of various values of uc relative to experimentally relevant 
timescales (e.g. j^^^^^^^ ^ different forms of the 

dynamics above cjc, and examine in detail the effects on 
qubit coherence. 



A. Spin Echo Decay 

Beginning with the spin echo we study the impact of 
uc and the form of the high-frequency rolloff on the mea- 
sured qubit coherence. Experimental data for coherence 
as a function of free-precession time have been shown to 
be fit well to a superexponential decay ex exp{ — {t/T2)^)^ 
with 7 = 4. We numerically calculate the spin-echo decay 
as a function of free-precession time and fit to a super- 
exponential for various values of uc^ and different high- 
frequency rolloffs. As seen in Fig. [2]3, varying ujc has a 
significant effect on the best-fit values of 7 extracted from 
the numerical calculations. As uc increases we observe a 
reduction of 7 = 4 ^ 3, corresponding to the presence of 
increased high-frequency spectral weight in Si3{uj). For 
the smallest values of uJc the 1/e coherence time ex- 
pected for a spin echo experiment increases dramatically 
beyond what has been observed. Offsetting these changes 
to the predicted j'j^^^^^^ by increasing a would be incon- 
sistent with microscopic modeling and would yield values 
of T2^^^^ inconsistent with experiment, thus limiting the 
range of possible values of uoc- These calculations suggest 
a possible value of cc;c/27r ^ 10-100 kHz to best approxi- 
mate experimental measurements. Such values are in line 
with analytical calculations of the high-frequency cutoff 
for nuclear-spin diffusion processes. 

The form of the high-frequency cutoff in Sjsiuo) also 
has significant impact on the experimentally measured 
qubit coherence. A microscopic, quantum mechanical 
treatment predicts that ODD and UDD sequences can 
deco uple the electron- nuclear interaction to arbitrary 
ordei^^^^. These results already indicate that (to the ex- 
tent that a description in terms of a noise spectral density 
is valid) the spectral density should have a high frequency 
cutoff that is faster than any power law. A uo~^ cutoff 
would lead to a exp(— t'^^^) decay for UDD or CDD de- 
coupling sequences of sufficiently high order according to 
a perturbative short-time expansion. We argue that a 
Gaussian rolloff where the noise power spectrum is sup- 
pressed as Sj^iuj) (X uj~'^e~^^^^^^ is most appropriate for 
the dynamics of nuclear spin diffusion, motivated by the 
the pair correlation approximation of Ref. 42 (see Ap- 
pendix for details). 

We present numerical simulations assuming this form 
of high-frequency rolloff, as well as power law scaling 
cj"^, corresponding to a rolloff of 3( dB/octave (e.g. 12 
dB/octave = uj~^). The transition from 7 = 4 ^ 3 as a 



function of uJc varies with the form of the high-frequency 
rolloff (Fig. [|d). Setting a;c/27r = 10 kHz in Fig. we 
show full numerical calculations for accumulated error 
in a spin-echo experiment as a function of free-precession 
time, superexponential fits to these data, and superexpo- 
nential fits with 7 fixed to be four, for different values of 
the high-frequency rolloff. We observe that for the steep- 
est values of the rolloff the calculated error rates are fit 
well using 7^4, while reducing the steepness of the 
rolloff (corresponding to increasing the weight of high- 
frequency components in S^{uj)) suppresses the best-fit 
7 towards three. Calculations using a Gaussian rolloff 
show good agreement with a superexponential decay us- 
ing 7 = 4. As shown in panel b, however, these differ- 
ences are reduced further as uJc 100 kHz. 



B. Multipulse Dynamical Decoupling 

Divergences between the superexponential fit with 7 = 
4 and full numerical calculations show differences that 
may be observable in experiments, but differences in the 
form of the spin-echo decay are subtle and will likely be 
limited by experimental noise. The scaling of qubit co- 
herence times and error rates in the high-fidelity regime 
with pulse number, n, provide further insight into the 
high-frequency nuclear dynamics. 

In Fig. [3] we show the calculated error rates in col- 
orscale as a function of free-precession time and n, for 
CPMG and UDD. We vary the value of uc and the form 
of the high-frequency rolloff between the "hard" Gaussian 
and the "soft" 12 dB/octave power-law to serve as perfor- 
mance bounds. In the high-fidelity regime constant-error 
contours demonstrate the efficiency of error suppression 
for times short compared to T2. As loc is reduced calcu- 
lated coherence times increases in all circumstances. In 
the case of a hard Gaussian rolloff, the form of the UDD 
filter function^ ^ provides significant benefits relative to 
CPMG and provides much deeper error suppression as 
well as enhanced coherence times. By contrast, the pres- 
ence of a slow 12 dB/octave rolloff suppresses any benefits 
of using UDD relative to CPMG. 

Aside from the absolute magnitudes of the extracted 
coherence times, the high-frequency components of S 13(00) 
also impact the form of the scaling of the calculated error 
suppression with n. For small values of uc the addition 
of pulses produces an approximately linear increase in co- 
herence time with n, as high-frequency fluctuations con- 
tribute little dephasing. For uc ^ I — 10 kHz, coherence 
times in the ms regime are possible for n ^ 20. For larger 
values of uJc coherence times of this order are possible 
with increasing n, but as discussed previously, the form 
of Sf3{uj) yields diminishing returns with n. The highest 
frequency components of interest are set approximately 
by r~l^^ the minimum interpulse free-precession time. 
When r~l^ ujc/^tt^ reduction in Tmin (via increasing 
n) fails to provide additional measurable suppression of 
noise. However, in the case of relatively large values of 
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uc small improvements in measured error rates may be 
obtained with increasing n. 

These characteristics are captured in the power-law 
scaling of Tf oc n'^'""""' (Tf ^^"^ oc n^'""""') 

Fig. [3j)-p shows the effect of the high-frequency cutoff 
on ip. For both CPMG and UDD lowering the value of 
ujc increases ^ 1. However, this phenomenon is most 
pronounced in the case of a hard Gaussian rolloff; the 
presence of a soft high-frequency rolloff suppresses the 
changes in with ujc 



C. Low- Frequency Nuclear Dynamics and 
Programming 

For the sake of a transparent and self-contained pre- 
sentation, we have assumed that the spectrum is well 
described by a power law down to a suitable low- 

frequency cutoff. A full treatment of nuclear spin diffu- 
sion shows that this approximation breaks down below 
'^0.1 Hz; below this frequency the power spectral density 
shows significant deviations from l/cj^ scaling. However, 
these frequencies only affect the relation between a and 
rj^{FiD) ^ whereas the behavior at higher frequencies rel- 
evant for decoupling sequences remains unaffected. The 
direct determination of a as described above provides 
an independent estimate that is not affected by the spe- 
cific form of the low-frequency spectral density. Further, 
calculations including modified low-frequency behavior 
show that the best-fit value of a remains in the range 

From the calculations presented heretofore it is readily 
apparent that one of the most effective ways to improve 
qubit coherence is to mitigate the fluctuating overhauser 
field experienced by the qubit. Nuclear spin diffusion has 
been characterized and even controlled^^^ through feed- 
back mechanisms, allowing the time-ensemble- averaged 
T2* to be extended. The most important parameter in 
such an experiment is the bandwidth of the noise sup- 
pression procedure; mitigating only low- frequency noise 
will not improve coherence limits with n-pulse dynami- 
cal decoupling sequences where the limiting performance 
is governed by rapid fluctuations at frequencies near the 
dynamical decoupling limit. For Sf3{(jj) oc l/cj^ we find 
that using experimentally relevant parameters, suppres- 
sion of the low-frequency tail has a dominant effect ex- 
clusively on T2^^^^ unless the bandwidth of the noise 
suppression extends beyond 100 kHz. As such, program- 
ming nuclear dynamics to mitigate low-frequency fluc- 
tuations can make T2^^^^ ji{Echo) ^ -g^ contrast, in 
the presence of a nuclear spectrum with a sharp cutoff, 
if the bandwidth of nuclear stabilization far exceeds the 
cutoff frequency, Si3{uj) is suppressed over its entire rel- 
evant range, as given by the product of Sf3{uj) and the 
filter transfer function defined by the nuclear program- 
ming protocol. In this case, high-bandwidth nuclear sta- 
bilization has the potential to dramatically extend coher- 



ence and improve the efficacy of dynamical decoupling 
sequence performance. 

V. NONINSTANTANEOUS CONTROL PULSES 

Realistic experiments are conducted with control tt 
pulses of nonzero duration (r^^ 7^ 0), mitigating the util- 
ity of the so-called bang-bang limit in which most dy- 
namical decoupling studies are performed. In previous 
work, a simple modification of the filter function was 
implemented in order to account, to lowest order, for 
the effect of noninstantaneous control pulsed. In this 
procedure, we modify the time-domain filter function by 
incorporating a delay with length r^^ and value zero be- 
tween free-precession periods. This approximation as- 
sumes that dephasing is negligible during the application 
of a TT pulse, consistent with many experimental observa- 
tions. Incorporating this delay results in a filter function 
for an arbitrary n-pulse sequence 

F{ujr) = \yn{uJr)\^ 

n 

= |1 + (-i)^+ie^^^ + 25] [-lye'^o^r (.Qg (a;r^/2)|2 

(7) 

where 8jT is the time of the center of the j*^ tt pulse, and 
r is the sum of the total free-precession time and 7r-pulse 
times. In previous experiments we showed that this ap- 
proximation could account for evolution of the noise field 
during the control pulses, improving the accuracy of the 
theoretical model. 

We have employed this model in the current setting to 
understand the effect of the finite bandwidth of the ap- 
plied pulses in experiments on the singlet /triplet qubit. 
Modeling r^^ G [0.1,20] ns, we find negligible effects on 
qubit coherence (changes < 10% in the calculated 1/e 
coherence time) for all spectra of interest (not shown). 
These changes are largest for high values of the high- 
frequency cutoff {uoc ^ 100 kHz), and typically less than 
1% for ujc ~ 10 — 30 kHz. Further, the finite bandwidth 
of the pulses did not produce discernible changes in the 
relative performance between the CPMG and UDD se- 
quences. Typical deviations in the high-fidelity regime 
demonstrated that the minimum accumulated error could 
be impacted at the 10~^ level for r^^ ~ 100 ns and 
uJc > 50kHz. These differences are not detectable in 
current experiments, suggesting that at present the ef- 
fects of pulse-bandwidth are negligible. 

The model presented above assumes ideal, but non- 
instantaneous control pulses. A significant source of er- 
ror in dynamical decoupling experiments may be derived 
from imprecise control operations. The presence of a 
large error even for small n in experiments on S/T qubits 
suggests that pulse-fidelities may be quite poor. One ex- 
pects that the effects of control pulse imperfection can be 
elucidated by experimental studies of error scaling with 
n; short-time behavior would show large error increases 
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FIG. 3: The effect of a high-frequency cutoff in the l/cj^ noise power spectrum on qubit coherence under the appUcation 
of multipulse dynamical decouphng. Qubit error probabihty (coherence) is expressed as a colorscale, with a value of 0.5 
corresponding to complete dephasing. Left panels correspond to CPMG, right panels correspond to UDD. Individual columns 
for panels (c-n) denoted with high-frequency rolloff in use. Black dashed line corresponds to the 1/e decay contour. White dashed 
lines correspond to constant error contours beginning at 1% and diminishing by 100 x per step, from right to left. Contours 
thus correspond to error probability 10~^, 10~^, 10~^, etc. Contours for values of error probability below 10~^^ are not shown, 
(o-p) Extracted T2 for CPMG (o) and UDD (p) given different values of ujc and the high-frequency rolloff. Lines are power-law 
fits to the data holding the n = intercept fixed at tI^^^^ — 26 ns, and varying V^. Extracted fit parameters: Gaussian Rolloff 
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these calculations, a = 5 x 10 
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with n in the presence of substantial pulse errors that 
increase with each applied operation. Examining the rel- 
ative performance of CP vs CPMG multipulse spin echo 
could also serve to illuminate the role of imperfect control 
pulses. 



VI. DISCUSSION AND CONCLUSION 

In this manuscript we have presented a phenomeno- 
logical model for the error accumulation in FID, Spin 
Echo, and multipulse dynamical decoupling experiments 
on Singlet-Triplet spin qubits in GaAs. Our simula- 
tions are based on an experimentally validated model 
for error accumulation in the presence of a noisy envi- 
ronment characterized by an arbitrary Si3{uj). The re- 
sults presented herein suggest that a simple semiclassical 
noise model can accurately reproduce a variety of experi- 
mental measurements, and can provide strong agreement 
with more detailed theoretical calculations. Further, they 
have elucidated some of the performance limits one might 
expect based on relevant noise processes associated with 
nuclear spin dynamics. 

The salient characteristics of experiments^^ providing 



the best measured values of Tr 



(Echo) 



and T: 



{CPMGr^) 



are 

reproduced by simulations incorporating S^{uj) cx I/cj^ 
in our phenomenological model. An extraction of of the 
relevant noise strength, a based on comparison with these 
measurements and T2^^^^ agrees within a factor of or- 
der unity with more detailed calculations of nuclear-spin 
diffusion. Further, our calculations provide data sup- 
porting the presence of a high-frequency cutoff ujc ^ 100 
kHz with Gaussian rolloff, consistent with heuristic mi- 
croscopic theory and experimental measurements. 

The results presented here provide predictive and an- 
alytical power for experimental studies of multipulse dy- 
namical decoupling, permitting detailed characterization 
of the form of S^{uj) via observation of the scaling of 



n{CPMGr^) 



(X 



j(CPMG) 



in 



(UDDr,) 



jAUDD) 



). Calcu- 



lations incorporating a high-frequency cutoff uJc ~ 100 
kHz are consistent with experimental measurements^^ 
of -0 ~ 0.7. Detailed information on the form of the 
high-frequency rolloff can be obtained by comparison of 
the performance of CPMG and UDD in the high-fidelity 
regime. Unfortunately, at this time measurement and 
operational infidelities mask the effects of decoherence- 
induced error accumulation to be studied in this regime. 

Achieving the ultimate benefits associated with UDD 
or other optimized dynamical decoupling protocols will 
thus require significant improvements in operational fi- 
delity, measurement fidelity, and classical noise filtration. 
The dynamics of nuclear spins likely provide the domi- 
nant dephasing mechanism in Singlet-Triplet qubits re- 
alized today, but as experimental capabilities improve, 
other noise sources will come into play. If, for instance, 
it is demonstrated clearly that nuclear spin dynamics ex- 
hibit a high-frequency cutoff with a Gaussian rolloff, for 



frequencies > cjc, an additional noise floor will likely 
come into play due to other processes (e.g. ambient mag- 
netic field fluctuations due to current noise near the sam- 
ple). Simulations show that including such a noise floor 
is effectively the same as modifying Sj3{(jj) to include a 
soft rolloff above the value of uoc set by nuclear spin dy- 
namics. These high-frequency spectral components limit 
coherence and reduce the benefits associated with the 
use of optimized decoupling strategies such as UDD. Un- 
derstanding system performance with soft noise rolloff 
therefore serves as a simple approximation for the ulti- 
mate influence of residual extrinsic noise processes. We 
have also shown that detailed, systematic studies of S/T 
qubit coherence under the application of various pulse 
sequences and values of n can provide a useful diagnostic 
for relevant noise sources, as an indirect form of noise 
spectroscop^^^. 

These simulations have demonstrated that the most 
effective way to improve qubit coherence is through sup- 
pression of high-frequency noise. The difference in cal- 
culated coherence times and error rates associated with 
soft and hard cutoff frequencies motivates effort in hard- 
ware engineering in order to suppress any effects of exter- 
nal noise sources. For instance, in the presence of weak, 
high- frequency-dominated electrical noise, improved sup- 
pression arising from a transition between single-pole and 
multi-pole filters in an experimental system could lead 
to coherence-time extension by 2-3 x in a dynamical de- 
coupling sequence. Limiting dephasing noise to intrinsic 
sources could more importantly suppress error rates by 
orders of magnitude in the high-fidelity regime via use of 
UDD relative to CPMG. 

Historically, a large focus of the community has been 
on the realization of spin qubits in materials dominated 
by zero-nuclear-spin isotopes such as carbon or silicon. 
The studies presented here support this general view- 
point as we believe the dynamics of nuclear spins in the 
GaAs host material to be the dominant source of ob- 
served dephasing. However, alternate noise sources such 
as fluctuations in the effective exchange coupling due to 
electrical voltage noise may in fact dominate such ex- 
periments, and will occur irrespective of the presence 
of a nuclear spin bath. If access to high-fidelity single- 
qubit operations is readily available, dynamical decou- 
pling pulse-sequence application may form an effective 
means to reduce performance gaps between different ma- 
terials systems, extending coherence by orders of magni- 
tude relative to T2^^^\ 



Appendix: Gaussian High-Frequency Rolloff 

The dynamics of nuclear spins are driven by the Knight 
shift due to the hyperfine interaction with the electron 
spin and interactions with neighboring nuclei. The latter 
drive flip-flop (i.e. spin transfer) processes between nu- 
clei of the same species, which over long time scales cause 
the spin diffusion discussed above. Inter-species flip-flops 
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are suppressed by the mismatch in Zeeman energies due 
to the apphed field. We argue that transitions of a sin- 
gle nuclear spin are driven by the transverse component 
of the (effective) interaction fields generated by nearby 
nuclei of the same species. The distribution of this field 
should be similar to the longitudinal component reflected 
in the NMR resonance line, but a factor 2-3 smaller since 
only the same-species resonant contributions are relevant 
here. Based on the measured NMR line shape and the 
fact that many randomly oriented spins contribute, we 
assume a Gaussian distribution. 

The Knight shift due to the hyperfine interaction with 
the electron spin is on the order of 10 G. However, its ef- 
fect on the dynamics is much smaller since the electronic 
wave function extends over many unit cells so that nearby 
nuclei experience approximately the same Knight shift. 
Only the Knight shift difference is relevant because it 
detunes the flip-flop interaction. Taking a typical length 
scale of 5 nm over which the wave function changes ap- 
preciably along the z-direction, the Knight shift variation 
between nuclei that are one lattice constant apart is on 
the order of 1 G, comparable to the transverse field set- 
ting the off diagonal matrix elements. Neglecting it is 
thus reasonable for the purpose of developing a rough 
phenomenological model. 

Furthermore, it should be noted that the change of the 



Overhauser field from each such a spin fiip-fiop is much 
smaller than that resulting from flipping a single nuclear 
spin, since the angular momentum is only transfered over 
short distances to other nuclei with approximately the 
same hyperfine coupling constant. Instead of using the 
above estimate of this reduction based on the wave func- 
tion shape, we simply approximate the Gaussian distribu- 
tion of flip-flop frequencies acts as a multiplicative factor 
on the diffusion spectrum discussed above. This phe- 
nomenological approach automatically solves the prob- 
lem that our argument is not valid at longer time scales, 
when the diffusive process associated with many subse- 
quent flip-flops sets in. At the corresponding low fre- 
quencies, the Gaussian factor is near unity, so that the 
diffusion spectrum is recovered. 
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